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We investigate domain wall junctions in a generalized Wess-Zumino model with a Zjv symmetry. 
We present a method to identify the junctions that are potentially BPS saturated. We then use a 
numerical simulation to show that those junctions indeed saturate the BPS bound for N — 4. In 
addition, we study the decay of unstable non-BPS junctions. 



I. INTRODUCTION. 

In field theories with discrete, degenerate vacua, do- 
main walls - field configurations that interpolate between 
different vacua - may occur. Supersymmetric theories, 
for which the vacuum energy is positive semi-definite, 
frequently have degenerate vacua. In many instances, 
multiple vacua exist in which the vacuum energy vanishes 
and supersymmctry remains unbroken. The degenerate 
vacua may or may not be related by a spontaneously 
broken discrete symmetry. 

General considerations concerning domain walls in su- 
persymmetric theories are presented in jl]-^. A lower 
bound on the tension, the BPS bound, can be calculated 
without explicitly solving for the profile of the wall. Do- 
main walls break the translation symmetry in one direc- 
tion, and supersymmetry is either completely broken, or 
1/2 of it is preserved. In the latter central exten- 

sion appears in the supersymmetry algebra. The tension 
saturates the BPS bound and is equal to the (1, 0) central 
charge. 

In it was noted that more complicated field config- 
urations with axial geometry, domain wall junctions, may 
also occur. These configurations, of the "hub and spoke" 
type, are a natural generalization of domain walls, and in- 
terpolate between more than two degenerate vacua. Far 
away from the center of the junction, the fields are ap- 
proximately constant - at the vacuum expectation values 
of the various vacua - in sectors. These sectors are sep- 
arated by "spokes" , where the field interpolates between 
one vacuum to the next following approximate domain 
wall profiles. In the center of the junction, "the hub" , the 
domain walls meet and the field configuration resembles 
a string. If these junctions are BPS saturated, 1/4 of the 
original supersymmetry is preserved and both the (1,0) 



and (1/2,1/2) central charges appear in the supersym- 
metry algebra. The junction tension is determined by a 
combination of the (1/2,1/2) and (1,0) central charges 
H. Although the two central charges contain (related) 
ambiguities, it is shown in that the ambiguities can- 
cel in physical quantities like the tension. 

It is worthwhile to emphasize that domain walls in su- 
persymmetric models are not necessarily BPS saturated. 
For example, in the model of || the BPS saturation of 
a domain wall depends on the values of the model pa- 
rameters. In many models it is not possible to obtain 
explicit analytical solutions to the first order BPS equa- 
tions. In those models the question whether solutions 
exist at all can still be answered. In a necessary and 
sufficient condition which indicates whether two vacua 
are connected by a PBS saturated domain wall was de- 
veloped^. 

The issue of BPS saturation in junctions is more com- 
plicated, because a similar condition does not exist. At 
this time, there is only one known model [jlOj in which 
a domain wall junction was found as an explicit ana- 
lytical solution to the BPS equations. In the absence 
of analytical expression for the field configuration, al- 
ternative methods must be employed to study junctions. 
Some information can be obtained from analytical results 
for specific limiting values of parameters or coordinates. 
Complementary results can be obtained from numerical 
simulation of the equations of motion. 

For any domain wall junction, a BPS bound on the 
junction tension can be calculated. In M the general 



*It is not obvious that this method works also for theories 
with more than one chiral superfield. 
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method to obtain this bound from the domain walls sur- 
rounding the junction is presented. In this work, it was 
also established that the junction tension is negative in 
general, although "exotic" models in which domain wall 
junctions have positive tension may exist. 

In this paper we will study domain wall junctions in a 
class of generalized Wess-Zumino (WZ) models with Zn 
symmetry. The Lagrangian is given by 

C = Jd 2 9 J d 2 6 + ^Jd 2 0W(<S>) + h.c. 

(1.1) 

where the Kahler potential has the canonical form 
/C(<I>, = , and the superpotential is given byf] 



W($) = $ 



1 



TV + 1 



N+l 



(1.2) 



The action is invariant under the transformation 
$(x, 0) -> e 2 ™/ N $(x, 9e- m / N ). Moreover, the model has 
TV physically equivalent vacua given by 



_ ^2-rrik/N 



fc = 0,l,...,TV- 1, 



(1.3) 



where <fi is the scalar component of the superfield $. De- 
spite its apparent simplicity, the model is of physical in- 
terest because it is related to the Veneziano-Yankielowicz 
action jl]J , which is a low energy effective action for su- 
persymmetric QCD. The parameter ./V is related to the 
number of colors in the SU(iV) gauge group. 

For TV = 2, the model reduces to the renormalizable 
Wess-Zumino model. It is well known that there is an an- 
alytical expression for the domain wall in this case (see 
M, for example). In the large TV limit, an analytical 
expression for the "basic" domain wall - the wall that 
connects vacua for consecutive values of k — is presented 
in |l^,[ll| to next to leading order in 1/TV. In this ref- 
erence it is also noted that all vacua in the model are 
connected by BPS saturated domain walls at finite TV. 

The "basic" domain wall junction has also been stud- 
ied. We define the "basic" junction as the field configura- 
tion that cycles through all TV vacua, consecutively from 
k = to k = N— 1, counter-clockwise around the center 
of the junction. In |lj] an analytical solution to the BPS 
equation for the "basic" junction in the TV — > oo limit 
was presented. In |j it was confirmed that the "basic" 
junction is also BPS saturated for finite values of TV by 



numerical simulations. The tension of the "basic" junc- 
tion was also calculated, analytically to next to leading 
order in 1/TV, and numerically also for finite values of N. 

Beside the "basic" junction, the vacua of the model 
allow for a multitude of other types of junctions. As dis- 
cussed further in the paper, we devise a general method 
to identify all potential BPS saturated junction in this 
model. There is a very stringent (asymptotic) consis- 
tency condition which allows only for a well defined set 
of intersection angles between the walls. The set of po- 
tential BPS junctions we identify contains the junctions 
that appear in 15 lq] , where tilings of domains and net- 
works of domain walls in this model are studied. 

We consider four types of domain wall junctions in the 
model with N — 4, two of which belong to the class 
of potential BPS saturated, and two more that are not 
BPS saturated. For each junction we perform a numer- 
ical simulaton of the second order equations of motion 
to obtain the junction profile. For the potential BPS 
junctions, we compare the energy of the junction to the 
BPS bound (which we evaluated from the domain walls 
surrounding the junction by the methods in 0) to de- 
termine whether they are indeed BPS saturated. For the 
other two junctions, we investigate their stability and de- 
cay. Although we perform our numerical calculations for 
the value TV = 4, we expect our results to apply to other 
values as well. Moreover, the methods we use are quite 
general, and can be applied to other models. 

T he re mainder of the paper is organized as follows. In 
Sec. [A, we describ e the pos sible domain walls of the 
model. Then in Sec. I A 1 and 1 A 2, we describe in detail 



the two types of domain walls which exist in the case 
N — 4. These walls appear far a away from the center of 
the junctions that we consider. In Sec. [ B we first discuss 
the potential BPS junctions in the models for generic N, 
which satisfy a very restrictive consistency condition. We 
then focus on the TV = 4 case and describe in detail the 



two types of BPS jun ctions (Sec. IB 1 and IB 2) and the 



non-BPS types (Sec. [B 3 ) which occur. In Sec. IC 



we 



give a brief outline of the computational scheme used for 
obtaining the numerical representation of the junctions; 
the paper ends with some comments in Sec. 



A. Domain walls. 

BPS saturated domain walls are solutions to the equa- 
tion 



t More generally we could consider the superpotential 
W($) = A$ - B$ N+1 /(N + 1), but for our purposes this 



equivalent to Eq. (1.2) by rescaling the field and the 
j. . i a i d\T75v j. „_j _ . A(i-m/N d-i/at 



In 



dinates as <j> -> (A/By 7 *^ and x„ -» a (1 ~ n)/n B~ 1/N 
the literature, the large TV limit of the model is studied with 
A = TV and B = 1/N N ~ 1 . 



(1.4) 



with the boundary condition that (j) approaches vacuum 
values at x — > — oo and x — > oo. Such solutions exist for 
particular values of the phase 5, which are determined by 
the equation 

l^{e' lS W{<j){x -> -oo))} = Im{c- lS W{(l)(x -► +oo))} . 

(1.5) 
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For the model with N = 4, there are two types of 
domain walls. The first type connects adjacent vacua, 
for which k differs by one. The second type connects 
opposite vacua, for which k differs by two. We will discuss 
each of these types in turn. 

1. Type I domain walls. 

Although it is not possible to find an analytical so- 
lution for the domain wall profile, the characteristic fea- 
tures can be readily uncovered. For definiteness, we focus 
on the wall that co nnec ts the vacua labeled by k = and 
k = 1. From Eq. ( |1.5[ ) it follows that for this wall the 
value of 5 is either — 7r/4 or 3ir/4. For 6 = —tt/A (3tt/4) 
the wall connects the vacuum with k = 1 (k = 0) at 
x = — oo to the vacuum with k = (k = 1) at x = oo. 

A general solution to the BPS equation contains one 
integration constant, indicating the freedom to shift the 
center of the wall. We fix this constant by choosing the 
center of the wall to be at the origin: The wall is then 
invariant under the transformation <p(x) — i<f>(—x). To 
study the wall close to the center, it is useful to write the 
field in terms of polar variables, <f> — pe la . The symmetry 
of the wall dictates that a — tt/4 and d x p — at the 
center. It then follows from the BPS equation that, for 
5 = -tt/4, 

d x a\ x=Q = -1±A, (1.6) 
pa 

where po = p(0). The value of po can be calculated by 
evaluating the constant of motion along the wall, i.e. 
Im {e~ lS W ((/))} at x — oo and at x = 0. As a conse- 
quence, po is the sole real root of the equation 

p + f = ^/2. (1.7) 

Numerically, po — 0.55514. The tension of the wall can 
be calculated exactly, and depends only on the value of 
(j) at the two ends of the wall; it is given by 

Q 

T[ I] = 2 \W((j>(x -> oo)) - W(4>(x -> -oo))| = -V2. 

(1.8) 

In Table I we present all domain walls of type I with their 
corresponding values for the phase <5. 
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TABLE I. Type I walls with their corresponding values of S. 

2. Type II domain walls. 

For definiteness, we focus on the wall connecting the 
vacua labeled by k = 2 and k — 0. For this wall, the value 
of 5 is either or it. For 6 = (6 = %) the wall connects 
the vacuum with k — 2 (fc = 0) at x — > — oo to the 
vacuum with fc = 0(fc = 2)ata;— >oo. The wall is purely 
real. Again, we fix the integration constant so that the 
wall is centered at the origin; with this choice, the wall 
is invariant under the transformation <f>(x) = —<f>(—x). 

The field vanishes at the origin, its derivative being 
(for 6 = 0) 

M a =o = !. (1-9) 
while the wall's tension in this case is 

t[ I] = 2 \W((f>{x -» oo)) - W{<p(x -» -oo))| = -|. 



(1.10) 



Table II lists all domain walls of this type with their 
corresponding value of 5. 
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3tt/2 
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TABLE II. Type II walls with their corresponding values 
of 8. 



B. Junctions. 

BPS junctions are solutions to the equation 

3^ = ^(1-^). (1.11) 

Here z = x + ty and d z = l/2(d x — id y ). Junctions 
are static configurations that have the geometry of a hub 
and spokes system. In the sectors between the spokes, 
the field 4> approaches the vacuum values; far away from 
the hub, <j) follows domain wall profiles from one sector 
to the next on trajectories perpendicular to the spokes. 
Each spoke is therefore associated with a domain wall. 
The hub is the center of the junction where the domain 
walls meet. 

A few considerations are important to consider what 
junctions are possible, and which of those junctions are 
BPS saturated. 

First, the sum of the forces exerted on the junction by 
the domain walls that are attached to it has to vanish for 
a static configuration. 
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Second, the value of the phase S in the BPS equation 
for the junction merely determines the orientation of the 
junction. The effective value of the phase for the BPS 
domain wall equation for a wall associated with a spoke 
that makes an angle a with the y axis is equal to 5 — a. 
Once the direction of one spoke is fixed, for example to 
point in the y direction, the value of S is fixed. For the 
remainder of the spokes, only specific values of a are then 
possible. This restriction on the possible values of a turns 
out to be a very stringent consistency condition. In the 
case of the model with N = 4, inspection of Tables I 
and II shows that this condition allows only two types 
of junctions as potential BPS candidates. We call these 
two types A and B respectively, and we investigate them 
below. 

Once it has been established that a static junction ex- 
ists, another interesting issue to determine is whether or 
not it is stable. If a junction is BPS saturated, it is guar- 
anteed to be stable (also see the discussion in 0). For 
non-BPS junctions, the issue is not so clear; if the bound- 
ary conditions also allow a BPS saturated configuration, 
then the non-BPS junction is at most meta-stable. 

One can actually devise a general method for generat- 
ing all the potential BPS saturated junctions in our class 
of models. In fact, given the presence of the N degener- 
ate vacua (1.3), one has that for N even (odd) there are 
N/2 ((AT - l)/2) types of domain walls. These different 
types of walls can be labeled by the integer I = \ki — k 2 \ 
or I = N — \k% — &2|j whichever is smaller (here k% and 
&2 indicate the vacua that are connected by the wall); 
moreover the tension of a wall, labeled by I, is equal to 



T (i) 



2\ / 2nl 

\ 2 - 2cos . 

N+ 1 V N 



(1.12) 



Now, we can represent each type of wall by a vector 
with length equal to its tension; thus a junction will be 
represented by the vectors of the walls that it contains 
far away from the origin, with the direction of the vectors 
chosen so that they point in the direction where the cor- 
responding wall is located. In a sense, the vectors play 
the role of the spokes in the hub and spoke system, but 
they have additional information because of their length. 
Crossing a spoke counter clockwise means that the value 
of k increases by I from one domain to the next along 
the corresponding domain wall: Thus, two junctions will 
be equivalent if they can be mapped onto each other by 
inversions and rotations. 

All of the potential BPS saturated junctions can then 
be obtained as follows. We start with the basic Ad- 
junction, which consists of N vectors representing walls 
with I = 1 (type I walls) at angles of 2-k/N . Then we add 
any two neighboring vectors in the diagram representing 
the junction, repeating this process until all uncquivalcnt 
diagrams have been generated. For AT = 4 this procedure 
gives the type A and B junctions of above (Fig. |l|). In 
Fig. [| all the potential BPS junctions are shown for the 
case N = 6. 




FIG. 1. Potential BPS saturated junctions for the case 
TV = 4. 




A 




FIG. 2. Potential BPS saturated junctions for the case 
N = 6. 



1. Type A junction. 

Here we study the junction of the type shown in 
Fig. |l| a). Far away from the center of the junction, there 
are four sectors in which the field consecutively takes on 
the four possible vacuum expectation values. These sec- 
tors are separated by BPS domain walls of type I, which 
were discussed before. 

Near the center of the junction, the field can be ex- 
panded in a power series in z and z. Accordingly one 
has 



(z,z) 



n — rn— 



(1.13) 



with a mn some complex coefficients 



'One can in principle allow also for negative n and m with 
the condition m > \n\ (n > \m\) if n < (m < 0). These 
negative va lues a re however excluded when o ne sub stitute the 
expansion (1.13) in the differential equation (1.11). 
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The differential equation (1.11) and the boundary con- 
ditions are invariant under the transformations (7 = 
-Tr/2 + vr/iV) 



Z 2 : 
Zjv : 



(z,z) 
(z,z) 



0t( e - 



211 z e 2 * 7 



e 2i7 z,-e- 217 ^). 



(1.14) 
(1.15) 



Imposing the Zjv symmetry the expansion ( 1.13| ) takes 
the form 



oc 00 



fc=0 1=0 



(1.16) 



with bki and c^z again complex coefficients. Imposing 
in addition the Z2 symmetry, we find the following con- 
straints on these coefficients 



b U = (-ire 2 ^ ;i 



Ckl 



(-1) 



JV(fe+l)„2j7_t 



(1.17) 
(1.18) 



Thus, if N is even, one can introduce the real coeffi- 
cients dki and eu, obtaining the expansion 



it. In order to determine whether those asymptotic solu- 
tions can be connected consistently and thus to investi- 
gate whether the junction is BPS saturated we performed 
a numerical simulation of the second order equations of 
motions for the field 0, while imposing the domain walls 
as the boundary conditions at large distance from the ori- 
gin. In order to allow the field to relax to its minimum 
energy configuration, we also introduced a damping term 
(see Sec. IC for details). In Fig. ||we show the modulus 
of <f> after the configuration has come to rest. 



. 5' 




2 . 5 



(Z,Z 



EE [dklZ l+Nk (zz) 1 +e kl z N{ - k+1) -\zz) 

(1.19) 



fe=0 1=0 



If on the other hand N is odd, one has 



b kl = {-l) k e^b\ 



kl- 



Ckl = (-i)( fc+ v* 7 cL 



kl' 



(1.20) 
(1.21) 



so that, according to the value of the index fc, half of the 
coefficients is real, the other half being purely imaginary. 
In the case N = 4 one has 7 = — 7r/4, so that 



(z,z) 



^/^\-[d kl z^ ik (z-z) l + e kl z wk {zz) 1 



k= o 1=0 



(1.22) 



We can now calculate the coefficients d k i and e k i using 
the BPS equation; it turns out that all coefficients are 
determined except for the ones which multiply terms that 
contain only z, that is, the e k o- These coefficients have to 
be adjusted to obtain the correct domain walls far from 
the center of the junction^. 

Close to the origin, the configuration looks like a string. 
We have proved then that there are sensible solutions to 
the BPS equations both near the origin and far away from 



§Note that the expansion up to terms of eleventh order in z 
and z contains only two adjustable parameters. 



FIG. 3. Modulus of <j> for the type A junction. 

In order to determine whether the junction is BPS sat- 
urated, we will first calculate the BPS bound on the en- 
ergy inside a square of "radius" R, with the center of 
the junction located at the center of the square. Assum- 
ing the junction to be BPS saturated the energy inside a 
square with radius R can be calculated entirely in terms 
of quantities related to the domain walls surrounding the 
junction when the radius R is much larger than the size 
of the junction. Explicitly, for large i?, the bound on the 
energy takes the form 



^BPS =T {A) +AT iI) K 

L A 1 



(1.23) 



Here L is a distance in the directi on pe rpendicular to 
the x-y plane. The second term in (1.23) is the contri- 
bution to the energy from the four type I domain walls 
connected to the junction; the first term represents in- 
stead the energy of the junction, which is equal to 



T, 



{A) _ 



a ■ dl, 



square 



where 



ai = Im {4>d Xi 4>} 



(1.24) 



(1.25) 



Now, each wall on the contour of the square contributes 
the same amount to the tension. The BPS bound on the 
junction tension is therefore equal to 



5 



T (A) = 4 



i {I) dx, 



(1.26) 



where w- 1 ' — Im dx^ 1 ^} an d <j>^ is the profile of a 
type I wall along the x axis. Using a numerical repre- 
sentation of the type I domain wall, we determined that 



T, 



U) 



-2.8614. 



We also calculated the energy in a square with radius 
R as a function of R directly from our numerical repre- 
sentation of the junction using the equation 



E 
L 



{d x $d x (t> + d y $d y <f> + V}dxdy, (1.27) 



square 



where V = d^Wd^W is the scalar potential. In Fig. || we 
compare this energy with the BPS bound. It is clear that 
the energy converges to the bound for large R (compared 
to the size of the junction) and therefore the junction 
does indeed saturate the BPS bound. 
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FIG. 4. Energy inside a square of radius R as a function 
of R (dots) compared to the BPS bound (solid line) for the 
type A junction. 



2. Type B junction. 

The second type of junction that is potentially BPS 
saturated has two spokes corresponding to type I walls at 
an angle of tt/2. A third spoke associated with a type II 
domain wall intersects each of the other spokes at an 
angle of 37r/4. The configuration is sketched in Fig. |l| b). 

This type of junction may be obtained from a junction 
of type A by pinching two neighboring spokes together 
in a symmetric fashion, thereby eliminating one domain. 
We performed a numerical simulation of the second order 
equations of motion, with the appropriate domain wall 
profiles as boundary conditions, and including a damping 
term (see again Sec. [IC] for details). In Fig. | we show the 
modulus of the field <j) after the configuration has come 
to rest. 



. 5' 
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FIG. 5. Modulus of 4> f° r the type B junction. 

To analyze whether this type of junction is BPS satu- 
rated we follow a procedure similar to the one used for 
the type A junction. We first calculate the bound on the 
energy in a square of radius R for large R, assuming the 
junction to be BPS saturated. This bound takes the form 



E 



BPS 



L 



= T, 



(B) 



iT^R - 



V2T{ n 'R — AT. (1.28) 



The first term represents the BPS bound on the tension 
of the junction. The second and the third term represent 
the contribution to the energy from the two type I walls 
and the type II wall (the factor of y/2 appears because 
the spoke representing the type II wall is in the diagonal 
direction). Special care must be taken because the spoke 
representing the type II wall intersects the sides of the 
square diagonally, so that the energy contributed by two 
triangular areas needs to be subtracted. This is the origin 
of the fourth term, AT, which takes the form 



AT = 4 / d x ^d x (S) {u ^xdx, 



(1.29) 



where (j)" 1 ' is the profile of a type II wall along the x 
axis. 

The field <f> on the type II wall connecting the k = 
sector to the k = 2 sector of this junction is purely real, 
so that it does not contribute to the BPS bound on the 
tension of the junction. On the other hand the tension 
receives equal contributions from each of the type I walls. 
The BPS bound on the tension of type B junctions is 
therefore half of the tension of a type A junction 



T, 



(B) 



1. 



-T, 



(A) 



(1.30) 



We can finally evaluate AT using the numerical repre- 
sentation of the type II wall, finding in this way 



AT 



1. 



-T, 



(A) 



(1.31) 



This is a curious equality, as AT is a quantity calculated 

(A) 

from a type II wall, and T 2 can be calculated from a 
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type I wall profile. We have no further comments on the 
nature of the equality here, but as a consequence the total 
energy in a square with large radius is equal for type A 
and type B junctions. 

We also calculated the energy in a square with radius 
R as a function of R directly from our nume rical rep- 
resentation of the type B junction using Eq. (1.27). In 
Fig. ^ we compare this energy with the BPS bound. It 
is clear that for this junction too, the energy in a square 
converges to the bound for large R (compared to the size 
of the junction). The type B junctions therefore also 
saturate the BPS bound. 




FIG. 6. Energy inside a square of radius R as a function 
of R (dots) compared to the BPS bound (solid line) for the 
type B junction. 



3. Non-BPS junctions. 

Here we discuss static junctions that are not BPS satu- 
rated. For such junctions the values of S associated with 
the various spokes are not consistent. However, the to- 
tal force on the junction vanishes. Moreover, they are 
invariant under an extended symmetry group, and there- 
fore extremize the energy. 

The first junction of this type that we discuss is dis- 
played in Fig. a). It consists out of four sectors sepa- 
rated by spokes which intersect at ir/2 angles. However, 
in contrast to the BPS four-junction, the spokes are as- 
sociated with type II walls. 



. 5' 
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FIG. 8. Modulus of 4> f°r the non-BPS four-junction. 

However, this junction has a negative mode. It is un- 
stable against local perturbations that break its symme- 
try, and it decays into a configuration with three domains 
separated by two domain walls, as shown in Fig. ^. Two 
domains with the same vacuum expectation value con- 
nect, and the resulting two domain walls are pushed out 
of the center. 
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FIG. 9. Modulus of <f> after the non-BPS four-junction has 
decayed into a configuration with three domains separated by 
two walls. 



I""--., 




FIG. 7. Non-BPS saturated junctions in the case N = 4. 

We found by a numerical simulation of the second order 
equations of motion that a static configuration of this 
type does indeed exist. In Fig. || we show the modulus 
of 4> for the final configuration at rest. 



The second type of non-BPS walls we discuss is the 
eight-junction shown schematically in Fig. [?] b). Com- 
pared to the BPS four-junction, this junction has winding 
number two. We again found by numerical simulation of 
the second order equations of motion that a static con- 
figuration of this type exists. We show the modulus of 
the field <j> for this configuration in Fig. [l(| 
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FIG. 10. Modulus of <f) for the non-BPS junction with wind- 
ing number two. 

We also observed that this junction can decay into two 
four-junctions with winding number one when a sym- 
metry breaking perturbation is introduced. The two re- 
sulting junctions repel each other, as shown in Fig. |ll| 
We noticed that the eight-junction only decays when the 
perturbation is sufficiently large. For now, we leave open 
the question whether this indicates that the junction is 
meta-stable, or that this behavior is an artifact of the 
finite size of the lattice we used. 




FIG. 11. Modulus of <j> after the non-BPS eight-junction 
has decayed into a configuration with two type A junctions. 

C. Numerical methods. 

We simulated the second order equations of motion 
on a lattice using a forward predicting algorithm. The 
lattice spacing was chosen to be much smaller than the 
size of the junctions, and, at the same time, the size of 
the lattice was much larger than the size of the junctions. 
For our purposes, a lattice of 251 by 251 points offered 
sufficient resolution. 

We added a damping term to the equations of motion 
which allows the field to relax to a minimum energy con- 



figuration. At the same time, this damping stabilizes the 
forward predicting scheme. 

In order to generate unstable junctions, we started our 
simulation with an initial configuration that was invari- 
ant under the same symmetry transformations as the 
boundary condition. As these symmetry are not broken 
by the equations of motion, the final configuration is the 
lowest energy configuration consistent with the symme- 
tries. Such a configuration may be unstable against local 
perturbations which break the symmetry. 

II. SUMMARY 

In this paper we have studied in detail domain wall 
junctions in a generalized Wess-Zumino model. We have 
presented a method to identify all potentially BPS sat- 
urated junctions, and we have described a procedure to 
determine whether these junctions indeed satisfy the BPS 
bound. We showed that in the case N = 4 (the lowest 
value of N for which there is more than one potential BPS 
junctions), these junctions are in fact BPS saturated. 

On the basis of our results, in conjunction with the 
results for the basic junction for generic value of N fij 
and the large N results for the basic junction B, we 
speculate that every potential BPS saturated junction for 
any N in this model indeed saturates the BPS bound. 
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